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ABSTRACT 



Aims. Based on a rather complete model atom for neutral and singly-ionized iron, we evaluate non-local thermodynamical equilibrium 
(non-LTE) line formation for the two ions of iron and check the ionization equilibrium between Fe i and Fe n in model atmospheres 
of the cool reference stars. 

Methods. A comprehensive model atom for Fe with more than 3 000 measured and predicted energy levels is presented. As a test and 
first application of the improved model atom, iron abundances are determined for the Sun and five stars with well determined stellar 
parameters and high-quality observed spectra. The efficiency of inelastic collisions with hydrogen atoms in the statistical equilibrium 
of iron is estimated empirically from inspection of their different influence on the Fei and Fen lines in the selected stars. 
Results. Non-LTE leads to systematically depleted total absorption in the Fei lines and to positive abundance corrections in agreement 
with the previous studies, however, the magnitude of such corrections is smaller compared to the earlier results. Non-LTE corrections 
do not exceed 0.1 dex for the solar metallicity and mildly metal-deficient stars, and they vary within 0.21 dex and 0.35 dex in the 
very metal-poor stars HD 84937 and HD 122563, respectively, depending on the assumed efficiency of collisions with hydrogen 
atoms. Based on the analysis of the Fe i/Fe n ionization equilibrium in these two stars, we recommend to apply the Drawin formalism 
in non-LTE studies of Fe with a scaling factor of 0.1. For the Fen lines, non-LTE corrections do not exceed 0.01 dex in absolute 
value. This study reveals two problems. The first one is that gf— values available for the Fei and Fen lines are not accurate enough 
to pursue high-accuracy absolute stellar abundance determinations. For the Sun, the mean non-LTE abundance obtained from 54 Fe i 
lines is 7.56±0.09 and the mean abundance from 18 Fen lines varies between 7.41+0. 1 1 and 7.56±0.05 depending on the source of 
the gf— values. The second problem is that lines of Fe i give, on average, a 0. 1 dex lower abundance compared to those of Fe n lines 
for HD 61421 and HD 102870, even when applying a differential line-by-line analysis relative to the Sun. A disparity between neutral 
atoms and first ions points to problems of stellar atmosphere modelling or/and effective temperature determination. 

Key words. Atomic data - Atomic processes - Line: formation - Stars: atmospheres - Stars: fundamental parameters 



1. Introduction 

Iron plays an outstanding role in studies of cool stars thanks to 
the many lines in the visible spectrum, which are easy to de- 
tect even in very metal-poor stars. Iron serves as a reference 
element for all astronomical research related to stellar nucle- 
osynthesis and the chemical evolution of the Galaxy. Iron lines 
are used to derive basic stellar parameters, i.e. the effective 
temperature, T e g, from the excitation equilibrium of Fei and 
the surface gravity, log g, from the ionization equilibrium be- 
tween Fe i and Fe n. In stellar atmospheres with 7^ > 4500 K, 
neutral iron is a minority species, and its statistical equilib- 
rium (SE) can easily deviate from thermodynamic equilibrium 
due to deviations of the mean intensity of ionizing radiation 
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from the Planck function. Therefore, since the beginning of 
the 1970s, a large number of studies attacked the problem of 
non-local thermodynamic equilibrium (non-LTE) line forma- 
tion for iron in the atmospheres of the Sun and cool stars. The 
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(120051) . and they were w idely applied in stellar parameter and 
abundance analyses (see lAsplundll2005l for references). It was 
understood that the main non-LTE mechanism for Fe i is ultra- 
violet (UV) overionization of the levels with excitation energy 
of 1.4 to 4.5 eV. This results in an underpopulation of neutral 
iron where all Fe i lines are weaker than their LTE strengths, and 
it leads to positive non-LTE abundance corrections. 

The need for a new analysis was motivated by the following 
problems uncovered by the previous non-LTE calculations for 
iron. 
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First, the results obtained for the populations of high- 
excitation levels of Fe i were not always convincing. The highest 
levels presented in the model atom did not couple thermally to 
the ground state of Fe n indicating substantial term incomplete- 
ness. To force the levels near the continuum into LTE, an up- 
per level thermalization procedure was applied (more d etails in 
IGehren et al.ll2001aHKorn e7alll2003t ICollet et al.ll2005l) . 

A second aspect is the treatment of poorly known inelas- 
tic collisions with hydrogen atoms. Their role in establish- 
ing the statistical equi librium of ato m s in cool stars is de - 
bated for decades, from lGehrenl d 19751) to lBarklem et al.1 d2010l) . 
Experimental data on H i collision cross-secti ons are only avail- 
able for the reson ance transition in Na i dFleck et al.l 1 1 99 it 
iBelvaev et al.l l 1999b. and detailed quantum mechanical calcula- 
tions were published for the transitions between first n ine levels 
in Li i (IBelvaev & Barklemll2003HBarklem et al.ll2003h and Na i 
( IBelvaev et al.1 1 19991 l2010t iBarklem et al.1 12010) . For all other 
chemical species, the basic formula us ed to ca l culate collisi ons 
with H i ato ms is the one proposed by iDrawinl d 19681 1 19691) . as 
described bv lSteenbock & Holwegerld 19841) . and it suggests that 
their influence is comparable to electron impacts. The labora- 
tory measurements and quantum mechanical calculations indi- 
cate that the Drawin formula overestimates rate coefficients for 
optically allowed transitions by one to seven orders of magni- 
tude. Therefore, various approaches were employed in the liter- 
ature to constrain empirically the efficiency of H i collisions. The 
studies of stellar Na i lin es favor a low efficiency of this type of 
collisions. For example, lAllende Prieto et ail d2004l) found that 
the center-to-limb variation of the solar Na 1 6160A line is repro- 
duced in the non-LT E calculations with pure electron collisions. 
iGratton etalJd 19991) calibrated H i collisions with sodium using 
RR Lyr variables and concluded that the Drawin rates should 
be decreased by two orders of magnitude. Based on t heir so- 
lar Na i line profile analysis, IGehren et alj d2004l) and iTakedal 
d 19951) recommended to scale the Drawin rates by a factor Sh 
= 0.05 and 0.1, respectively. With a similar value of Sh = 0.1, 
the ionization equilibrium between Cai and Can in selected 
metal-poor stars was matched consist ently with surface gravi - 
ties derived from Hipparcos parallaxes dMashonkina et al.l2007l) . 
On the other hand, spectroscopic studies of different chemical 
species suggested that the H i collision rates might be reason- 
ably well described by Drawi n's formula wi t h Sh > 1 ■ For ex 



ample, empirical estimates by IGratton et aL ( 19991) resulted in 
Sh — 3 for O i and Mg i and Sh = 30 for Fe i. Allende Prieto et a"fl 
d2004l) and lPereira etail d2009l) inferred Sh = 1 from the analysis 
of the center-to-limb variation of solar O i t riplet ~ 7770A lines. 
The same value of Sh = 1 was obtained bv ITakedal d 19951) from 
solar Oi line profile fits. For a review of stu dies cons t rainin g 
empirically the e fficie ncy of H i collisions , see iLambertl d!993l) : 
iHolwegerid 19961) . and lMashonkinal d2009l) . 

As a result of applying incomplete model atoms and a differ- 
ent treatment of collisions with hydrogen atoms, no consensus 
on the expected magnitude of the non-LTE effects was achieved 
in the previous studies of iron, and res ults were in confli ct with 
each other in some cases. For example. iKorn et ail d2003l) found 
a negligible discrepancy between the non-LTE spectroscopic and 
Hipparcos astrometric dista nces of the halo star HP 84 937, while 
the non-LTE calculations of lThevenin & Idiarnd 1999h resulted in 
a 34 % smaller spectroscopic distance of that same star. 

This study aims to construct a fairly complete model atom 
of iron, to be tested using the Sun and selected cool stars 
with high-quality observed spectra and reliable stellar param- 
eters. Compared with the previous non-LTE analyses of iron, 



the model atom of Fe i was extended to high-l ying levels pre- 
dicted by the atomic structure calculations of iKuruczl d2009l) 
and this turned out to be crucial for a correct treatment of the 
SE of iron in cool star atmospheres. With our improved model 
atom we tried to constrain the scaling factor Sh empirically. We 
realize that the real temperature dependence of hydrogen col- 
lision rates could be very different from that of the classical 
Drawin formalism, and we may not always achieve consistent 
Sh values from the analysis of different stars. We also realize 
that the required thermalizing process not involving electrons in 
the atmospheres of cool metal-poor stars could be very differ- 
ent from inelastic collis i ons w ith neutral hydrogen atoms. For 
example. IBarklem et al._d2pi0|) (see also Barklem et al.l 120031 : 
IBelvaev & Barkleml 120031 : Belvaev et al.l 1201 oh uncovered the 
importance of the ion-pair production and mutual neutralisation 
process A(nl) + H( 1 s) ^ A + + H for the SE of Li and Na. Since 
no accurate calculations of either inelastic collisions of iron with 
neutral hydrogen atoms or other type processes are available, 
we simulate an additional source of thermalization in the atmo- 
spheres of cool stars by parametrized H i collisions. Investigating 
the Sun as a reference star for further stellar differential line-by- 
line analysis, we also derive the solar iron absolute abundance 
and check the solar Fe i/Fe n ionization equilibrium using an ex- 
tended list of lines, which can be detected at solar metallicity 
down to [Fe/H] = -2.5. We find it important to inspect the accu- 
racy of atomic data for various subsamples of iron lines in view 
of comprehensive abundance studies across the Galaxy targeting 
at stars of very different metallicities. 

This paper is organized as follows. The model atom of iron 
and the adopted atomic data are presented in Sect. [2] There we 
also discuss how including the bulk of predicted Fe i levels in 
the model atom affects the SE of iron. In Sect.fJJ the solar iron 
spectrum is studied to provide the basis for further differential 
analyses of stellar spectra. Section|4] describes observations and 
stellar parameters of our sample of stars, and Sect. [5] investi- 
gates which line-formation assumptions lead to consistent ele- 
ment abundances from both ions, Fe i and Fe n. Uncertainties in 
the iron non-LTE abundances are estimated in Sect.[6] Our rec- 
ommendations and conclusions are given in Sect.|7] 



2. The method of non-LTE calculations for iron 

In this section, we describe the model atom of iron and the pro- 
grams used for computing the level populations and spectral line 
profiles. 



2.1. Model atom 

Energy levels. Iron is almost completely ionized throughout the 
atmosphere of stars with an effective temperature above 4500 K. 
For example, nowhere in the solar atmosphere, the fraction of 
Fei exceeds 10%. Such minority species are particularly sen- 
sitive to non-LTE effects because any small change in the ion- 
ization rates changes their populations by a large amount. To 
provide close collisional coupling of Fe i to the continuum elec- 
tron reservoir and consequently establish a realistic ionization 
balance between the neutral and singly-ionized species, the en- 
ergy separation of the highest levels in the model atom from the 
ionization limit must be smaller than the mean kinetic energy of 
electrons, i.e., 0.5 eV fo r atmospheres of sola r temperature. 

In our earlier study (IGehren et al.ll2001al hereafter Paper I), 
the model atom of Fe i was built up using all the ener gy lev- 
els from the experimental analysis of iNave et al. I (fl99l . in to- 
tal 846 levels with an excitation energy, E sxc , up to 7.5 eV. The 
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Fig. 1. The atomic structure of Fei as obtained from laboratory experiments (long lines, black lines in the color online version of 
the paper) and calculations (short lines, blue lines in the color online version of the paper). See text for the sources of data. 



lat er updates of ISchoenfeld et al.l (1 19951) and the measurements 
of IBrown et al. I (11988b provided another 112 energy levels for 
Fe i. All the known levels of Fe i are shown in Fig.[T]by long hor- 
izontal bars. How complete is this system of l evels? In the sa me 
figure, the upper energy levels predicted by iKuruc 3 (120091) in 
his calculations of the Fei atomic s tructu re are plotted by short 
bars. As described by iGrupp et"ai1 d2009l) . the new calculations 
for iron included additional labo ratory levels an d more config- 
urations than the earlier work of lKuruczl (119921) . It can be seen 
from Fig. Q] that the system of measured levels is complete be- 
low £exc - 45 000 cirT 1 (5.6 eV), except perhaps for singlets. 
However, laboratory experiments do not see most of the high- 
excitation levels with £ exc > 7.1 eV, which should contribute a 
lot to provide close collisional coupling of Fe i levels to the Fe n 
ground state. 

In our present study, the model atom of Fe i is constructed us- 
ing not only all the known energy levels, but also the predicted 
levels with £ exc up to 63 697 cirT 1 (7.897 eV), in total 2970 
levels. The measured levels belong to 233 terms. Neglecting 
their multiplet fine structure we obtain 233 levels in the model 
atom. The predicted and measured levels with common parity 
and close energies are combined whenever the energy separation 
is smaller than 150 cirT 1 at E exc < 60000 cirT 1 and smaller than 
210 cirT 1 at E exc > 60000 cirT 1 . The remaining predicted lev- 
els, all above E exc = 60000 cirT 1 , are used to make up six super- 
levels. For super-levels, the energy is calculated as a ^-weighted 
mean and the total statistical weight amounts from 940 to 2160. 
Our final model atom of Fe i is shown in Fig. [2] 

For Fe n (Fig.|2]i, we rely on the reference model atom treated 
in Paper I and use the levels belonging to 89 terms with £ exc up 
to 10 eV. Multiplet fine structure is neglected. The ground state 
of Fe m completes the system of levels in the model atom. 

Radiative bound-bound (b-b) transitions. In total, 11958 
allowed transitions occur in our final model atom of Fei. 
Their aver a ge-"m ultiplet" ^/-values are calcu lated us i ng the 
iNave et al.l J 1994 compilation for 2649 lines, iKuruczl d2009l) 



calculations for the transiti ons between th e measured levels, in 
total, for 73 434 lines, and IKuruczl J2009I) calculations for the 
transitions between the measured and predicted and between two 
predicted levels, in t otal, for 2 81 00 7 lines. The quality of the re- 
cent cal culations of Kuruczl d2009l) is estimated by comparing 
with the lNave et alJ ( 1994 ) data. The latter can be referred to as 
ex perimental data th ough, for a minority of the lines compiled 
by INave et al.l d 19941) . ^/-values were derived from solar spec- 
tra assuming a solar iron abundance. Taking the experimental 
sample as reference the theoretical data show a single line scat- 
ter of 0.54 dex. Ignoring all lines with deviations above 1 dex 
(6 %) the scatter is reduced to 0.33 dex. The reliability of the 
calculated ^/-values is thus roughly characterized by a factor of 
2 statistical accuracy only. It is, however, important to note that 
there is no systematic shift between experimental and calculated 
data. The mean difference is A log ^/(calculated - experimental) 
= -0.04 for the whole sample of lines. For 1525 allowed transi- 
tions in Fe 11, ^ /-values are fully based on the data calculated by 
IKuruczl dl992h . 

Radiative bound-free (b-f) transitions. Photoionization is the 
most important process deciding whether the Fei atom tends 
to depart from LTE in the atmosphere of a cool star. As in 
Paper I, our non-LTE calculations r ely on the phot oionization 
cross-sections of the IRON project dBautistal Il997l) . They are 
available for all the levels of Fe 1 in our model atom with the 
ionization edge in the UV spectral range, in total, for 149 levels. 
For the remaining levels of Fe 1 and all the Fe 11 levels, the hy- 
drogenic approximation was used. We note that photoionization 
weakly affects the SE of Fe 11 because Fe m constitutes only an 
extremely small fraction of the total iron atoms. 

Collisional transitions. All levels in our model atom are cou- 
pled via collisional excitation and ionization by electrons and by 
neutral hydrogen atoms. For Fe 1, electron-impact excitations are 
not yet known with sufficient accuracy, and, as in Paper I, our 
calculations of collisional ra tes rely on theoretical approxima- 
tions. We use the formula of Ivan Regemorterl d 19621) for the al- 
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Fig. 2. The model atom of Fe i (bottom panel) and Fe n (top panel). The super-levels are indicated by long lines. The iron spectral 
lines used in abundance analysis arise from the transitions shown as continuous lines. 



lowed transitions and assume that the effective collision strength 
T = 1 for the forbidden transitions. The latter assumption can be 
verified for the single forbid den transition of Fe i , a 5 D - a 5 F, 
for which the calculations of iPelan & Berringtonl (1 19971) in the 
IRON Project lead to T = 0.98 at T = 4000 K. 



For the transitions between the Fen terms up to z D°, 
we employ the data fro m t he close-coupling calculat i ons o f 
IZhang&Pradhanl (119951) and lBautista & Pradhanl (Il996lll998h . 
The same formulas as for Fe i are applied for the remaining tran- 
sitions in Fe n. Ioni zation by electr onic collisions is calculated as 
in Paper I from the Seatonl ( 1962) classical path approximation 



with a mean Gaunt factor set equal to g — 0.1 for Fe i and to 0.2 
for Fe ii. 

For collisions with H i atom s, we employ the formula 
of ISteenbock & Holwegeri (11984 for allowed b - b and b - 
f transitions and, following iTakedal ( 1 19941) . a simple rela- 
tion between hydrogen and electron collisional rates, Ch - 
C e ylJjnJmH)NH IN e , for forbidden transitions. The efficiency of 
Hi collisions is treated as a free parameter in our attempt to 
achieve consistent element abundances derived from the two ion- 
ization stages, Fe i and Fe n, in the selected stars. For each object, 
the calculations were performed with a scaling factor of 5h = 0, 
0.1, 1, and 2. In Fig. [3] we compare the electron-impact exci- 
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Fig. 3. Electron (open circles) and hydrogen (filled circles) col- 
lision rates of representative transitions in Fe i (left column) and 
Fe ii (right column) as a function of the transition energy in the 
model 5780/4.44/0 at a depth point of log T5000 = -0.5, where 
T = 5583 K, n H = 0.922 ■ 10 17 , and n e = 0.137 ■ 10 14 (top 
row), and in the model 4600/1.60/ - 2.50 at a depth point of 
logT 50 oo = -0.5, where T = 4174 K, n H = 0.212 ■ 10 17 , and 
n e = 0.286 • 10 11 (bottom row). Everywhere, 5k = 0.1 has been 
assumed. 



tation rates with the corresponding H 1 collision rates for repre- 
sentative transitions in Fe 1 and Fe 11 in the line-forming layers 
(logT 5 ooo = -0.5) of the solar-metallicity (r e ff/log g/[Fe/HU = 
5780/4.44/0) and metal-poor (4600/1. 60/-2.50) models. At the 
selected depth point, the electron number density n e in the solar- 
metallicity model is by a factor of 500 higher than that in the 
metal-poor and cool atmosphere, while the difference in neu- 
tral hydrogen number density iih is much smaller, only a fac- 
tor of 4. As a consequence, the hydrogen collisions dominate 
the total collisional rate for most Fe 1 and Fe 11 transitions in the 
metal-poor model. For solar-metallicity models, the influence 
of inelastic Hi collisions on the SE of iron is expected to be 
weaker because C///C e > 1 is fulfilled only for the Fei tran- 
sitions with small energy separation, i.e., smaller than 0.5 eV. 
For Fe 11, the close-coup ling electron-impact excitation rates of 
IZhang&Pradhanl(ll995h and lBautista & Pradhanl (1 19961 1 19981) 
are higher compared to the hydrogen collision rates even for 
small transition energies. A similar relation between Ch and 
C e also holds outside log 75000 = -0.5 in the solar-metallicity 
model. In the metal-poor model, the weakest iron lines form in- 
side log T5000 = -0.5, where the role of hydrogen collisions is 
weakened due to a decreasing «///«<> ratio. 



2.2. Programs and model atmospheres 

We compute statistical equilibrium populations for Fe 1 and Fe 11 
while keeping the atmospheric structure fixed. This is justified 
by the following considerations. Despite the fact that iron is 
an important source of the UV continuous opacity in an atmo- 
sphere of close-to-solar metallicity (see Fig.|4] top panel), the 
variations in its excitation and ionization state between LTE and 
non-LTE are found to have, at the stellar parameters with which 
we are concerned, no significant effect on the emergent fluxes 
(Fig-El bottom panel). Therefore, only minor effects on atmo- 
spheric temperature and electron number density distribution are 
expected. 

The solution of the non-LTE problem with such a compre- 
hensive model atom as treated in this study is only possible, 
at present, with classical plane-parallel (one-dimensional, ID) 
model atmospheres. All calculations were performed with model 
atmospheres computed with the MAFAGS-OS code dGruppI 
12004 lOrupp et al.ll2009h . which is based on up-to-date contin- 
uous opacities and includes the effects of line-blanketing by 
means of opacity sampling. 

We used a revised version of the DETAIL program 
dButler & Giddingsl[T98l based on the accelerated lambda it- 
er ation (ALI) scheme following the efficient method described 
bv lRvbicki & Hummed (Il99ll Il992h in order to solve the cou- 
pled radiative transfer and statistical equilibrium equations. The 
opacity package of the DETAIL code has been updated by 
the inclusion of the quasi-m olecular Lyman a satellit es fol- 
lowing the implem entation by ICastelli & Kuruczl (1200 ll) of the 
lAllard etaLl ( 1998 ) theo ry and by the use of the Opacity Project 
(see ISeaton et al.l 19941 for a general review) photoionization 
cross-sections for the calculations ofb-f absorption of C 1, N 1, 
Mg 1, Si 1, Al 1, and Fe 1. In addition to the continuous background 
opacity, the line opacity introduced by H 1 and metal lines was 
taken into account by explicitly including it in solving t he radi- 
ation t ransfer. The metal line list was extr acted from th e Kuruczl 
(11994 compilation and VALD database dKupkaet all 1 19991) . It 
includes about 720 000 atomic and molecular lines between 500 
and 300 000 A. The Fe 1 and Fe 11 lines as well as b -f absorption 
processes were excluded from the background. 

All the b - b and b - / transitions of Fe 1 and Fe n were 
explicitly taken into account in the SE calculations. The 57 
strongest b - b transitions were treated using Voigt profiles and 
the remaining transitions using depth-dependent Doppler pro- 
files. Microturbulence was accounted for by inclusion of an ad- 
ditional term in the Doppler width. 

The departure coefficients, b, = n? [LTE /nJ ;rE , were the n used 
to co mpute synthetic line profiles via the SIU program (iReetzl 
Il99lh . Here, «^ LTE and «[' TE are the statistical equilibrium and 
thermal (Saha-Boltzmann) number densities, respectively. In 
this step of the calculations, Voigt profile functions were adopted 
and the same microturbulence value V m j c as in DETAIL was ap- 
plied. 



2.3. Effect of model atom completeness on the statistical 
equilibrium of iron 

The departure coefficients for selected levels of Fei and Fen 
in two model atmospheres 5777/4.44/0 and 4600/1. 60/-2. 50 
are presented in Fig. [5] (top row). Our calculations sup- 
port the results of t he previous non - LTE s tudies of iron 
(lAfhav & Litesl Il972t iBovarchuk et all 1 1985b ITakedal 11991 



1 In the classical notation, where [X/H] 
log(N x /N H ) Slm . 



logtffc/iVk), 



[ Gratton et al.l 119991: Ifhevenin & Idlarll 119991: 
l2001at IShchukina & Truiillo Buenol 1200 it IShchukina et al 
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Fig. 4. Top panel: emergent fluxes calculated using the model of T e g = 5777 K, log g = 4.44, and [M/H] = with i ron (log sp R = 7.50, 
non-LTE with Sr =0.1, continuous curve) and without iron (dotted curve) compared to the solar flux taken from lWoods et al.l d 19961 
, dashed curve). The observed flux below 1700 A is contributed by the solar chromosphere. Bottom panel: the ratio between the 
non-LTE (Sh = 0.1) and LTE emergent fluxes in the same model atmosphere with logspe = 7.50 from the calculations with the 
final (continuous curve) and reduced (dotted curve, see Sect. 12. 3t model atoms. For clearer illustration, the theoretical fluxes were 
convolved with a Gaussian profile of AAd = 2 mA. 



120051: ICollet et alj|2005l) in the common stellar parameter range 
in the following aspects. 

- The non-LTE mechanisms for Fe i and Fe n are independent 
of effective temperature, surface gravity, and metallicity. 

- All the levels of Fe i are underpopulated in the atmospheric 
layers above log T5000 = 0. The main non-LTE mechanism 
is the overionization caused by superthermal radiation of a 
non-local origin below the thresholds of the levels with exci- 
tation energy of 1.4 to 4.5 eV. 

- Fe n dominates the element number density over all atmo- 
spheric depths. Thus, no process seems to affect the Fen 
ground-state and low-excitation-level populations signifi- 
cantly, and they keep their thermodynamic equilibrium val- 
ues. 

In this study, progress was made in establishing close col- 
lisional coupling of the Fei levels near the continuum to the 
ground state of Fe 11, due to including in the model atom a bulk 
of the predicted high-excitation levels of Fe 1. There is no need 
anymore in the en forced upper level thermalization proced ure 
that was applied bylGehren et all d2001al) : lKorn et al.l OOOlh . or 
ICollet etal.l(l2005l) . 

To illustrate the difference between the present and previ- 
ous non-LTE studies, we performed test calculations with the 
atomic model, where the predicted levels of Fe 1 are removed. 
Hereafter, the latter model atom is referred to as the reduced 



model. It includes 233 terms and 10740 radiative b - b transi- 
tios in Fe 1 and exactly the same levels of Fe 11 as in our final 
model atom. The bottom row panels of Fig. [5] show the depar- 
ture coefficients computed with the reduced model atom. It is 
evident that (i) the Fei terms are more underpopulated in the 
line-formation layers compared to the populations obtained with 
our final model atom (Fig. [5] top row), (ii) the highest levels do 
not couple thermally to the ground state of Fe n, and (iii) the ma- 
jority of the Fe 1 terms are close together. The explanation is that 
the net ionization-recombination rate of the Fe 1 levels separated 
by no more than 0.4 eV from the continuum is much lower in the 
reduced model atom than in the final model atom, for example, 
by a factor of 1000 and 150 at logrsooo - _ 0.6 in the model 
atmospheres 5777/4.44/0 and 4600/1. 60/-2.50, respectively. 

We found that emergent fluxes calculated with the 
MAFAGS-OS solar model atmosphere and using o ur final model 
atom reproduce well the observed solar fluxes of IWoods et al.l 
(fT9%l) in the 1700 - 3000 A spectral range, where neutral iron is 
an important source of the continuous opacity (Fig.|4j top panel). 
The observed flux below 1700 A is contributed by the solar chro- 
mosphere. As demonstrated in the bottom panel of Fig.|4] the 
difference in emergent fluxes between LTE and non-LTE is neg- 
ligible when employing the final model with Sh =0.1, while the 
use of the reduced model atom leads to a weaker b - f opacity 
of neutral iron compared to the LTE value, such that, in the 1600 
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Fig. 5. Departure coefficients, b, for the levels of Fe i and Fe n as a function of log 7-5000 in the model atmospheres 5777/4.44/0 (left 
column) and 4600/1. 60/-2. 50 (right column) from the calculations with our final model atom of iron (top row) and the reduced 
model atom which ignores the predicted levels of Fe i. Every fifth of the first 60 levels of Fe i is shown. They are quoted in the 
bottom right part of each panel. All the remaining higher levels of Fe i are plotted by dotted curves. For Fe n, every fifth of the first 
60 levels is shown. They are quoted in the top right part of each panel. The letters z, y, x, ... are used to denote the odd parity terms 
and a, b, c, ... for the even parity terms. To distinguish Fe i and Fe n levels, note that b < 1 for each Fe i level outside log T5000 — in 
the solar metallicity model and log T5000 = 0.4 in the metal-poor one. An exception is the highest levels in the layers with log T5000 
around -2.5 in the model 4600/1 .60/-2. 50. For Fe 11, b > 1 inside log T5000 = _ 2 in both model atmospheres. Everywhere, Sh = 0. 1 . 



- 2900 A spectral range, the change in emergent fluxes is about 
20%. 

In Sect.|6l we evaluate the difference in solar and stellar non- 
LTE abundances when using the final and reduced model atoms. 



3. Analysis of the solar Fe 1 and Fe 11 lines 

In this section, we derive the solar Fe abundance from the Fe 1 
and Fe 11 lines. The Sun is used as a reference star for further 
stellar abundance analysis. Th e solar flux observat ions are taken 
from the Kitt Peak Solar Atlas dKurucz et alll984l) . The calcula- 
tions were p erformed with the M AFAGS-OS model atmosphere 
5777/4.44/0 dGrupp et alJl2009l) . A depth-independent microtur- 
bulence of 0.9 km s _1 was adopted. Everywhere in this study, the 
element abundance is determined from line profile fitting. The 
uncertainty in the fitting procedure was estimated to be less than 
0.02 dex for weak and moderate strength lines (see Fig.|6| and 
less than 0.03 dex for strong lines. Our synthetic flux profiles are 
convolved with a profile that combines a rotational broadening 
of 1.8kms _1 and broadening by macroturbulence with a radial- 
tangential profile. The most probable macroturbulence velocity 
V mac varied mainly between 2.6 and 3.3 km s _I for different lines 
of neutral ir on and betwe en 3.4 and 3.8 km s _1 for Fe 11 lines. For 
comparison, lGravl(fT977l) found solar macroturbulence velocities 
varying between 2.9 and 3.8 km s for a small sample of the so- 



lar Fe 1 lines. The V mac values obtained from the non-LTE (5h 
= 0.1) fits of individual solar iron lines are indicated in Table[3] 
(online material). 

3.1. Line selection and atomic data 

The investigated lines were selected according to the following 
criteria. 

- The lines should be almost free of visible/known blends in 
the Sun. 

- For each star, the spectrum should include lines of various 
strength to provide an abundance analysis of both close-to- 
solar metallicity and very metal-poor stars. 

- The list of the Fe 1 lines should cover as large as possible a 
range of excitation energies of the lower level to investigate 
the excitation equilibrium of neutral iron. 

The selected lines are listed in Table|5] (online material), along 
with the transition information and references to the adopted 
^/-values. Van der Waals broadening of the iron lines is ac- 
counted for using the most accura te data available from c alcula - 
tions of lAnstee & O'Maral (I1995I): iBarklem & O'Maral dl997l): 
IBarklem et al.l dl998l) . and IBarklem & Aspelund-Johanssonl 
d2005h . 

Despite the existence of many sources of ^/-values for 
neutral iron, there is no single source that provides data for 
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Table 1. Solar mean iron non-LTE (Sh = 0.1) abundances deter- 



various sources 


of gf- 


-values 




Species 


N* 


log E 


cr 


Source 


Fei 


54 


7.56 


0.09 


Table[5] 


Fei 


54 


7.56 


0.09 


VALD 


Fen 


18 


7.56 


0.05 


Tabled 


Fen 


18 


7.45 


0.07 


VALD 


Fen 


18 


7.47 


0.05 


MB09 


Fen 


16 


7.56 


0.06 


M83 


Fen 


8 


7.41 


0.11 


SSK04 



* the number of lines. 
Sources: VALD = Kupka et al. (1999) . 
MB09 = Melendez & Barbuv (2009) . 
M83 = Moitv (1983) . 
SSK04 = Schnabel et al. (2004). 



all th e selected Fei lin e s. We employ e xperimental ^/-values 
fromjBard et alj (|1991l): iBard & Kockl d 19941): iBlackwell etalj 



(fl97 



1982aMbh : iFuhr et al.l dl988h . and lO'BrianetalJ dl99ll) . 
The accurac y of available sets of ^/ -values for Fei was 
discussed by iGrevesse & Sauval d!999l) . and the influence of 
using various sets of da ta on the de ri ved so lar iron abun- 
dance was inspected by iGehren et alj d2001bl) . Both papers 
recommen ded to employ the ^/-values publi shed by the 
Hannover (Sard et al.l jl99lt IBard & Kockl 11994 and Oxford 
dBlackwell et al.l ll979Ul982allbl) groups. Unfortunately, these 
groups did not provide data for some Fei lines important 
for the stellar iro n spectrum analysis . Five s ets of oscilla - 
tor strengths from [Melendez & Barbuvl d2009h: iMoitvl d!983l) : 
iRaassen & Uvlingsl (119981): ISchnabeletalJ d2004 . and the 
VALD database dKupka et al.ll 19991) were inspected for Fe n. We 
found that the solar mean abundance derived from the Fe n lines 
varies between logs = 7.41 ±0.11 and logs = 7.56 + 0.05 
depending on the source of g/- values (Table [TJ. Hereafter, the 
statistical error is the dispersion in the single line measure- 
me nts about the mean : cr = VS(x - x,) 2 /( w - 1). As reported 
by Grevesse & Sauval dl999l) . the data of IRaassen & Uvlingsl 
d 1998b for the visible Fen lines were underestimated by 
0.11 dex, on average, compared to the gf— values obtained us- 
ing lifetime and branching fraction measurements. 

3.2. Abundance analysis 

For each line, the element abundance was determined under var- 
ious line-formation assumptions: non-LTE with Sh = 0, 0.1, 1, 
and 2, and LTE. The quality of the fits is illustrated in Fig. [6] 
for two Fei lines. The results are shown in Table|5] (online ma- 
terial) and Fig. [7] To present the ab solute abundances f r om th e 
Fen lines we use gf— values from IRaassen & Uvlingsl (1 19981 ). 
where available, and lMoitvl (1 1 983b . just because they provide the 
mean element abundance consistent with that from the lines of 
Fe i. Compared with their LTE strength Fe i lines become weaker 
and Fen lines stronger under non-LTE conditions. The effect is 
largest for non-LTE with pure electron collisions (Sh = 0): the 
non-LTE abundance correction, Anlte = log £nlte _ log £lte, 
ranges between +0.03 and +0.15 dex for various Fei lines and 
between 0.00 and -0.02 dex for the Fen lines. As expected, the 
departures from LTE are weaker in the calculations with H i col- 
lisions taken into account. For the Fei lines, Anlte does not 
exceed 0.09, 0.04, and 0.03 dex, when S H = 0.1, 1, or 2, re- 
spectively. For the Fe n lines, the non-LTE abundance correc- 
tion is smaller than 0.01 dex in absolute values, independent of 




g 0.6 - 



5434.2 5434.4 5434.6 

Wavelength (A) 



1.00 - 



0.95 - 



0.90 - 



0.85 - 




6607.9 



6608.0 

Wavelength (A) 



Fig. 6. The best fits (continuous curve) of the solar Fei 5434 
and 6608 A lines (bold dots in top and bottom panels, respec- 
tively) obtained from the non-LTE calculations with Sh = 0.1. 
Fitting parameters are given in Table[5](online material). For Fe i 
6608 A, dotted curves show the effect on the theoretical profile 
of a ±0.02 dex variation in the final Fe abundance. The macro- 
turbulence velocity was reduced by 0.3kms~' for Aloge = 
-0.02 dex and increased by 0.2kms _I for Aloge = +0.02 dex. 
The bottom panel illustrates that errors due to an abundance vari- 
ation in fitting the solar iron line profiles are small. 



the Sh value. Our calculations showed that the low-excitation 
lines of Fei are subject to stronger deviations from LTE than 
the higher excitation lines, in agreement with the previous re- 
sults. For example, with Sh = 0.1, Anlte ranges between +0.04 
and +0.09 dex for six lines with E exc < 1 eV, while Anlte = 
+0.04 dex is the upper limit for the remaining Fei lines. As a 
result, an excitation energy gradient of the non-LTE iron abun- 
dances was found to be smaller than that for the LTE case. 
For example, with Sh = 0.1, we obtainecQ logepei = 7.516 + 
0.0154x£ exc (eV). 

We see no significant abundance discrepancies between low- 
excitation and high-excitation lines. However, the abundance 
scatter is considerably too high for E exc > 4 eV, independent 
of the line-formation assumptions used. For example, the differ- 
ence between the highest (Fe i 5517 A) and lowest (Fe i 5367 A) 
absolute solar abundance is about 0.5 dex in all the cases. The 
abundance difference between Fe i 5662 A and Fei 536 7 A with 
gf- values from a common source dO' Brian et al 1 199 11) amounts 
to about 0.25 dex. There are essentially two explanations for this 
result: either the ^/-values are not reliable or a significant frac- 



2 The subscript in log epci indicates that the iron abundance was de- 
termined from the neutral species. Likewise, a subscript Fe II refers to 
singly ionized lines used in the abundance determination. 
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Fig. 7. Solar LTE (left column) and non-LTE (5h = 0.1, right column) Fe abundances from the Fei (open circles) and Fen (filled 
circles) lines plotted as a function of E exc (top row) and W bs (bottom row). The mean abundances derived from the Fei and Fen 
lines along with their standard deviations are quoted in the top row panels. The dotted line indicates the mean abundance derived 
from the Fe i lines and the shaded grey area its statistical error. In the bottom row panels, the dashed line shows the calculated trend 
of log e Fe i versus W bs- Everywhere, gf- values indicated in Table[5]were used. 



tion of the iron lines are contaminated by relatively strong (yet 
undetected) blends or both. Note that, with an adopted micro- 
turbulence velocity of V m i c = 0.9 kms -1 , the trend of logepei 
versus W„b s is completely removed in the non-LTE calculations. 
The Fei 5517 A line is an obvious outlier in the plots of Fig. [7] 
its abundance is 0.3 dex higher than the mean of the remaining 
Fei lines. The solar mean Fei abundances indicated in Table[3] 
were computed without Fe i 5517 A and Fe i 5367 A taken into 
account, and the corresponding statistical error is large, <x = 
0.09 dex. These results demonstrate how far spectroscopic meth- 
ods lead us towards accurate and reliable absolute stellar (so- 
lar) iron abundance results. Some value around <xi og£ =0.1 dex 
seems to denote the current precision limit. 

The difference between the solar mean Fe i abundances cal- 
culated with Su = and 1 amounts to +0.07 dex. The corre- 
sponding value for Fe n is -0.01 dex. Since the amplitude of the 
non-LTE effect is smaller than the combined statistical error of 
the obtained Fe i and Fe n abundances and smaller than the dif- 
ference in Fe n abundances caused by employing various sources 
of gf— values, we did not endeavor to use the absolute solar Fe 
abundances determined from two ionization stages to constrain 
a 5h value. 



3.3. Comparison with other studies 

During the last decade possibly the most advanced approach to 
solar abundance determinations is based on ab initio 3D, time- 
dependent, hydrodynamical model atmospheres. We compared 
our ID LTE and non-LTE solar Fe abundance determinations 
with the results of 3D calculations performed independently by 
three scientific group s. The first 3D solar abundance analysis 
for Fe was made by lAsplund et all d2000l , ANTS2000). For 
nine common Fe i lines, the difference between their 3D+LTE 
and our 1D+LTE abundances amounts to, on average, -0.10 ± 
0.03 dex. The difference is smaller for the four common Fen 
lines when using the same ^/-values, A(ANTS2000 - this 
study) = -0.04 + 0.01 dex. In their recent 3D determinations, 



lAsplund et all (120091. AGSS2009 ) applied gf- values of the Fe n 
lines from Melendez & Barbuvl (l2009h and ID non-LTE cor- 
rection of +0.03 dex for Fei and got logepei = 7.52±0.05 and 
log epeii = 7.50±0.04 (here, the error bars is the systematic error 
added in quadrature with the statistical error calculated as the 
weighted standard error of the mean). Our ID non-LTE abun- 
dances loggpei = 7.56+0.09 (Table[3] Su = 0.1) and loggpen 
= 7.47 +0.05 when using ^/-values from IMelendez & Barbuvl 
(I2009h (Table[U are consistent within the error bars with the cor- 
responding values of AGSS2009, however, the abundance differ- 
ence (Fe i - Fe n) = 0.09 dex is larger compared to the 0.02 dex in 
AGSS2009. This could be due to the difference between ID and 
3D modelling. Note that updating the method of 3D calculations 
in AGSS2009 compared to ANTS2000 resulted in the smaller 
abundance difference (ID - 3D) for Fe i and opposite sign of the 
3D effect for Fe n. 

The first non-LTE stu dy for neutral iron that went beyon d 
the ID analysis was by IShchukina & Truiillo Buena d200ll) . 
Ignoring inelastic collisions with Hi atoms in their 1.5D+non- 
LTE calculations, they found a 0.074±0.03 dex higher solar Fei 
based abundance compared to the LTE value. Our lD+non- 
LTE analysis with pure electronic collisions gives for the so- 
lar Fei lines a very similar mean non-LTE abundance cor- 
rection of 0.08 dex. However, the absolute abundance is a 
0.11 dex higher compared to logg p e i = 7.50+0.10 obtained by 
IShchukina & Truiillo Buenol(l2001l) . 

Our ID Fe ii based abundance (the non-LTE effects are negli- 
gible) is 0.04 dex lower comp ared to loggpeii = 7.512+0.062 ob- 
tained bv lCaffau et alj ( 20101) from their solar 3D analysis when 



using ^/-values from IMelendez & Barbuvl (120091) in both stud- 
ies. 

Thus, the difference between our solar iron ID abundances 
and the data from the literature based on the 3D calculations 
does not exceed 0.04 dex when using either the Fe i or Fe n lines. 
However, the 3D effect is of opposite sign for the two ions, and 
this may affect analysis of the solar Fe i/Fe n ionization equilib- 
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Table 2. Characteristics of observed spectra. 



Object 


V 1 


Telescope / 


Observing run, 


Spectral range 


R 


S/N 




(mag) 


spectrograph 


observer 


(A) 






HD 10700 (r Cet) 


3.50 


2.2-m / FOCES 


Oct. 2001, K. Fuhrmann 


4500 - 6700 


60 000 


> 200 






2.2-m / FOCES 


Oct. 2005, F. Grupp 


4120 - 6700 


40000 


>200 


HD 61421 (Procyon) 


0.34 


2.2-m / FOCES 


Feb. 2001, A. Korn 


4200 - 9000 


60000 


>200 


HD 84937 


8.28 


VLT2 / UVES 


ESO UVESPOP 


3300 - 9900 


80000 


>200 


HD 102870 03 Vir) 


3.61 


2.2-m / FOCES 


May 1997, K. Fuhrmann 


4500 - 6700 


60000 


> 200 


HD 122563 


6.20 


VLT2 / UVES 


ESO UVESPOP 


3300 - 9900 


80000 


>200 



1 Visual magnitude from the SIMBAD database. 



4. Stellar sample, observations, and stellar 
parameters 

Our current sample consists of five stars with the effective tem- 
perature and surface gravity determined from methods largely 
independent of the model atmosphere. Two of them, HD 61421 
(Procyon) and HD 102870 (J3 Vir), are close-to-solar metallic- 
ity stars. HD 10700 (r Cet) is a mildly metal-deficient star with 
[Fe/H] =i -0.5. The two remaining stars are very metal -poor 
(VMP) with [Fe/H] < -2. HD 84937 represents the hot end of 
the stars that evolve on time scales comparable to the Galaxy 
lifetime. HD 122563, in contrast, is a cool giant. The selected 
stars are listed in Table[2] 

The spectroscopic observations for three of our program 
stars were carried out with the fibre-fed echelle spectrograph 
FOCES at the 2.2m telescope of the Calar Alto Observatory 
during a number of observing runs, with a spectral resolv- 
ing power of R 60000 (HD 10700 was also observed at 
R =i 40000 to get a higher signal-to-noise (S/N) spectrum at 
A < 4500 A). Characteristics of the observed spectra are sum- 
marized in Table[2] The S /N ratio is 200 or higher in the spectral 
range A > 4500 A, but lower in the blue. All the investigated iron 
lines lie longwards of 4427 A, and their profile analysis profits 
from the high spectral resolution and high S /N. For HD 84937 
and HD 122563, we used high-quality observ ed spectra from the 
ESO UVESPOP survey (iBagnulo et alJl2003h . 

Procyon, /3 Vir, and r Cet are among the very few stars for 
which the whole set of fundamental stellar parameters except 
metallicity can be determined from (nearly) model-inde pendent 
methods. In the most recent study of lBruntt et alj d2010h . the an- 
gular diameters of these stars were measured using interfero- 
metry. The measured bolometric flux, combined with the angu- 
lar diameter, implied r eff = 5383+47, 6494+48, and 6012+64 K 
for r Cet, Procyon, and f3 Vir, respectively. A surface gravity 
of log g = 4.54+0.02, 3.98+0.02, and 4.13+0.02, respectively, 
was comp uted using the stel lar mass and linear radius from 
Table 2 of lBruntt et alj d2010h . Our adopted stellar parameters 
of these three stars (TableO were taken from different studies, 
ho wever, they turned o ut to be consistent within 1 cr with the data 
of lBruntt et aLI d2010t) . We em ployed Procyon's ste l lar pa rame- 
ters from a careful analysis of lAllende Prieto et al.1 (120021) . The 
interferometric measurements of the angular dia meter and the 
bolometric flux were used by iNorth et all (120091) to derive T e g 
and log g of (5 Vir. F or t Cet, r e ff and log g were derived by 
iDi Folco et al.1 (|2004|) based on direct interferometric measure- 
ments of the angular diameter and stellar evolution calculations. 

For H P 84937, the i nfrare d flux method (IRFM) t emper - 
atures of lAlonso et all (Il996l): iMelendez & Ramfrezl (120041): 
iGonzalez Hernandez & Bonifacio! d2009l) . and lCasagrande et al.l 
d2010l) . r„ ff = 6330 K. 6345 K. 6333 K, and 6408 K, respectively, 



are consistent within lcr. Our adopted value, T e ff = 6350 K, 
is the average of the four. To determine log g, we used an 
updated Hipparcos p arallax of nnip = 13.74+0.78 mas from 
Ivan Leeuwenl d2007l) and a mass of 0.8 M derived from the 
tracks of lVandenBerg etail d2000l) . 

For HD 122563, we adopted r eff (IRFM) = 4600 K recom- 
mended by iBarbuv et al. (12003 ) an d based on the IRFM de- 
termination by lAlonso et al.1 d!999l). Using the 2MASS pho- 
tometric svstemT lGonzalez Hernandez & Bonifa cio (2009J) ob- 
tained a 200 K higher IRFM temperature of this star. However, 
they noted that the 2MASS photometric accuracy is very low 
for bright giant star s. The gravity was c alculated with nm p - 
4.22+0.3 5 mas from van Leeuwenl d2007l) and a mass of 0.8 M Q 
following IBarbuv et al.l ( 2003 ). The gravity errors indicated in 
Table|3] for HD 84937 and HD 122563 reflect the uncertainty in 
the trigonometric parallax. 



5. Fe i versus Fe n in the reference stars 

In this section, we derive the Fe i and Fe n abundances in the se- 
lected stars under various line-formation assumptions, i.e., non- 
LTE with Sh = 0, 0. 1 , 1 , and 2 and LTE, and we investigate which 
of them leads to consistent element abundances from both ion- 
ization stages. 



5.1. Methodology 

To minimize the effect of the uncertainty in gf— values on the fi- 
nal results, we applied a line-by-line differential non-LTE and 
LTE approach, in the sense that stellar line abundances were 
compared with individual abundances of their Solar counter- 
parts. Our results for Fe abundances are based on line pro- 
file analysis. In order to compare with observations, computed 
synthetic profiles were convolved with a profile that com- 
bines instrumental broadening with a Gaussian profile, rota- 
tional broadening, and broadening by macroturbulence with a 
radial-tangential profile. Rotational broadening and broadening 
by macroturbulence were treated separately only for Procyon 
and (3 Vir, with V sin / = 2.6 km s -1 and 2.5kms~', respectively 
dFuhrmannlll998l) . For the remaining stars, their overall effect 
was treated as radial-tangential macroturbulence. The macro- 
turbulence value VmaK w as determined for each sta r either by 
| Fuhrmannl dl998l l2004h or in our previous studies dKorn et al.l 
12.003b iMashonkina et al.ll2.d08h from the analysis of an extended 
list of lines of various chemical species. For a given star, V mac 
was allowed to vary by ±0.3kms~' (lcr). 

For each star, the microturbulence velocity was determined 
from the requirement that the iron non-LTE abundance derived 
from Fei lines with 5h = 0.1 must not depend on the line 
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Object 


Teg, K 


log g 


Ref. 


Vmic' 


Ion 


N 2 




Iron abundances, [Fe/H] 


















T TC 


Sr = 


5 H = 0.1 




1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


Sun 3 


5777 


4.44 




0.9 


Fe i 


54 


7.53±o.09 


7.61±0.10 


7.56±0.09 


7.54±0.09 












Fen 


18 


7.56±o.05 


7.55±o.05 


7.56±o.05 


7.56±0.05 


HD 10700 


5377 


4.53 


DTK04 


0.8 


Fei 


39 


-0.49±o.02 


-0.44±0.03 


-0.49±o.02 


-0.49±o.03 


(t Cet) 










Fe ii 


13 


-0.52±o.04 


-0.52±0.04 


-0.52±o.04 


-0.52±o.04 


HD 61421 


6510±49 


3.96±o.02 


AP02 


1.8 


Fe i 


45 


-0.12±o.07 


-0.12±o.07 


-0.10±0.07 


-0.13±o.07 


(Procyon) 










Fen 


15 


-0.03±o.05 


-O.O3±0.05 


-0.04±o.05 


-0.03±o.05 


HD 84937 


6350±37 


4.09±0.05 


IRFM+Hip 


1.7 


Fei 


21 


-2.17±o.07 


-1.96±o.06 


-2.00±0.07 


-2.13±0.07 












Fe ii 


9 


-2.08±0.04 


-2.06±0.03 


-2.08±o.04 


-2.08±0.04 


HD 102870 


6060±49 


4.11 ±0.01 


NDR09 


1.2 


Fe i 


39 


0.11 ±0.03 


0.11 ±0.04 


0.13±o.04 


0.11 ±0.03 


03 Vir) 










Fen 


13 


0.22±o.04 


0.23±0.04 


0.22±0.04 


0.22±0.04 


HD 122563 


4600±61 


1.60±0.07 


IRFM+Hip 


1.95 


Fei 


35 


-2.77±o.n 


-2.42±0.09 


-2.61±0.09 


-2.74±0.10 












Fen 


15 


-2.56±o.07 


-2.45±o.i2 


-2.56±o.07 


-2.56±0.07 



1 Microturbulence velocity is given in kms" 1 . 

2 Number of lines used. 

3 For the Sun, we present the absolute iron abundances calculated using g/-values from Table|5] 
References to the adopted r cff and log g: 

DTK04 = Pi Folco et al. (2004) . AP02 = Allende Prieto et al. (20021 . NDR09 = North et al. (20091 . 
IRFM+Hip = IRFM temperature and HippARCos-parallax based log g. 



strength. We note that, even for the other line-formation scenar- 
ios, the slope of the [Fe/H]i - W b s plot is also largely removed 
with the resulting V m i c value. 

Figure[8] shows our final non-LTE (5h = 0.1) abundances 
from individual lines of two ionization stages in each of the stars 
investigated. The average Fe i and Fe n abundances obtained in 
LTE and in non-LTE with 5h = 0, 0.1, and 1 are presented in 
Table [3] Having in mind the shortcomings of ID atmospheric 
structure modelling, we tend to represent for stars of various 
metallicity a similar range of the line formation depths. The 
strongest iron line detected in HD 122563 has W b s = 123.8 mA, 
while some of the selected iron lines are very strong in Procyon, 
B Vir, and r Cet with W ohs up to 250, 305, and 496 mA, respec- 
tively. Therefore, we used in stellar abundance and microturbu- 
lence velocity analyses only lines with equivalent widths less 
than 125 mA. 

5.2. Non-LTE effects 

5.2.1. Absolute abundances from Fen lines 

In the stellar parameter range with which we are concerned, 
the non-LTE effects of the investigated Fe n lines are found to 
be negligible when H i collisions are taken into account in the 
SE calculations. For the stars with metal abundances close to 
solar or mildly deficient, non-LTE with pure electronic colli- 
sions leads to a strengthening of the line core for the strongest 
(W > 40 mA) Fe ii lines because the line source function (S i u 
b u lbi B v ) drops below the Planck function (B v ) in the uppermost 
atmospheric layers. Here, b u and bi are the departure coefficients 
of the upper and lower levels, respectively. The lower levels of 
the investigated Fe ii lines are all of even parity with E exc = 2.8 
- 3.9 eV (see Table|5] online material), and they keep their ther- 
modynamic equilibrium populations (b/ = 1) throughout the at- 
mosphere. The population of the odd parity quartet and sextet 
terms with E exc = 5 - 6 eV is controlled by the strong transitions 
to the ground state a 6 D and low-excitation states a 4 F, a 4 D, and 
a 4 P. The line center optical depth of these transitions drops be- 



low unity at log 7/5000 < ~2 (Fig. [5), resulting in photon losses in 
the corresponding lines and underpopulation of the upper levels 
(b u < 1). For Procyon, B Vir, and t Cet, the non-LTE abun- 
dance correction is negative and does not exceed 0.02 dex in 
absolute value. The departures from LTE are negligible for the 
weaker lines. Non-LTE with Sh = results in the opposite ef- 
fect for the VMP stars. This can be understood as follows. In 
the atmosphere with [Fe/H] = -2.5, the line wing optical depth 
drops below 1 far inside the atmosphere even for the strongest 
Fe n transitions arising from the ground and low-excitation states 
to the odd parity levels with E exc =5-6 eV. As a result, they 
are pumped by the ultraviolet J v - B v {T e ) excess radiation and 
produce enhanced excitation of the odd parity levels in the lay- 
ers with logTsooo between +0.3 and -1 (Fig. [5]l. Here, the term 
"pumped" (or pumping) is used following iBruls et al.l dl992l) . 
In the upper layers, logTsooo < — 1> where the line center op- 
tical depth of the UV transitions drops below 1, enhanced ex- 
citation changes with photon losses resulting in a decrease of 
the departure coefficients of the odd parity levels. The investi- 
gated Fe n lines are formed in the atmospheres of HD 84937 and 
HD 122563 inside logTsooo = -1.4, where, for every line, the 
source function exceeds the Planck function, resulting in weak- 
ening the line relative to its LTE strength. The average non-LTE 
abundance correction amounts to +0.01 dex and +0.10 dex, re- 
spectively. 

5.2.2. Absolute abundances from Fei lines 

In all cases, non-LTE leads to a weakening of the Fe 1 lines. In 
close-to-solar metallicity models, this is mainly due to overion- 
ization. For each line, the source function is quite similar to the 
Planck function because most levels behave similarly, as can be 
seen in Fig. [5] (left top panel). At the other end of the metallic- 
ity range, the energy levels become weakly coupled far inside 
the VMP atmospheres due to deficient collisions (Fig. [5] right 
top panel). At the depths where the Fei lines are formed, the 
upper levels are all depleted to a lesser extent relative to their 
LTE populations than are the lower levels. The lines are weaker 
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Fig. 8. Non-LTE (5h = 0. 1 ) differential abundances derived from the Fe i (open circles) and Fe n (filled circles) lines in the selected 
stars as a function of observed equivalent width (left column panels) and excitation energy of the lower level (right column panels). 
In each panel, the dashed line shows the mean iron abundance determined from the Fe i lines. See text for more details. 



relative to their LTE strengths not only due to the general overi- 
onization but also due to b u /b/ > 1 resulting in S i„ > B v and the 
depleted line absorption. The non-LTE effects are the strongest 
when H i collisions are neglected in the SE calculations. On aver- 
age, Anlte(Sh = 0) = +0.05 dex and +0.06 dex for Procyon and 
P Vir and increases with decreasing metallicity to +0.13 dex at 
[Fe/H] = -0.5 (t Cet), +0.26 dex at [Fe/H] = -2.1 (HD 84937), 
and +0.43 dex at [Fe/H] = -2.5 (HD 122563). With H i collisions 
taken into account, the effect on the iron abundance determina- 
tion is substantially weaker. For example, with Sh = 0.1, Anlte 
= 0.00, +0.04, 0.03, 0.19, and 0.19 dex for the same sequence of 
stars, respectively. 



5.2.3. Differential abundances from Fei lines 



For stellar differential abundances, the non-LTE correction can 
be introduced as [Anlte] = [Fe/H]NLTE - [Fe/H]LTE- Since the 
departures from LTE for the solar Fe i lines are small in the cal- 
culations with Sh + 0, the corresponding [Anlte] and Anlte val- 
ues turned out to be close together. For example, [Anlte] (Sh = 
0.1) = +0.02, 0.02, 0.00, 0.17, and 0.16 dex for Procyon, p Vir, 
t Cet, HD 84937, and HD 122563, respectively. In the S H = 
case, the non-LTE differential abundance correction is approxi- 
mately 0.08 dex lower than the corresponding Anlte value. 
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5.2.4. Fei excitation equilibrium 

For the two VMP stars of our sample, non-LTE leads to the 
shallower excitation gradient of the Fe i line abundances com- 
pared to that obtained under the LTE assumption. For example, 
withSn = 0.1 (the right column panels ofFig.[8]l, d[Fe/H]/dE exc 
= +0.013 dex/eV for HD 84937 and -0.030 dex/eV for 
HD 122563, while the corresponding LTE values are -0.019 and 
-0.054 dex/eV, respectively. However, the non-LTE gradient ob- 
tained for HD 122563 is still steep compared to other stars of the 
sample. Our test calculations show that it could be substantially 
reduced, down to d[Fe/H]/dE exc = -0.008 dex/eV, with an ef- 
fective temperature revised downward by 80 K (l<x). No signif- 
icant change in the excitation slope between non-LTE and LTE 
was found for Procyon, /3 Vir, and t Cet, because the departures 
from LTE for Fe i are small in these stars. The gradient is small, 
at the level of -0.005 to -0.007 dex/eV, for r Cet. For Procyon 
and p Vir, d[Fe/H]/dE exc 0.02 dex/eV. We caution against de- 
termining stellar effective temperatures of metal-poor stars from 
the Fe i LTE excitation equilibrium. 

5.3. Empirically constraining the efficiency of H i collisions 

As can be seen from Table [3] the behavior of the difference 
between Fei and Fen abundances, A(Fen - Fei) = [Fe/H]n - 
[Fe/H]i, within various line-formation scenarios is different for 
stars of different metal abundance. 

A disparity between the neutral and singly-ionized iron is 
found in Procyon and in ft Vir, independent of either LTE or 
non-LTE. The least imbalance achieved in non-LTE with Sh 
0. 1 amounts to 0.06 dex and 0.09 dex, respectively. Very sim- 
ilar results for F e i/Fe n in Procyon were obtained earlier by 
lKornetalJd2003b . For the same star, a similar prob lem was un- 
covered also for Cai/Can dMashonkina et al.1l2007l) . Such a dis- 
parity can be caused by the uncertainties either in atmosphere 
and line-formation modelling or in stellar p arameters. The three- 
dimensional LTE simulations for Procyon (lAllende Prieto et all 
120021) showed that weak lines (W A < 50 mA) of both Fe i and 
Fe ii are weakened compared to a classical ID analysis, such that 
the derived Fe abundance of this star increases by 0.05 dex and 
0.04 dex, respectively. Thus, LTE+3D modelling does not mas- 
ter Fe i/Fe n in Procyon. The solution of the non-LTE+3D prob- 
lem for iron is still a cha llenge for s t ellar s pectra modelling (see 
Sect. [6]). As discussed bv lKorn et alj d2003l) . to remove a discrep- 
ancy of 0.07 dex between Fen and Fei in Procyon, one would 
require either a gravity of log g = 3.81, many standard devia- 
tions away from the astrometric result, or a temperature of T e g = 
6600 K, almost 2cr higher than the temperature based on bolo- 
metric flux. Our test calculations with an 80 K upward revised 
J e (f achieved consistent Fe i and Fe n abundances of this object 
(Tablegl. 

We find that non-LTE with Sh > 0.1 is as good as LTE to 
achieve the ionization equilibrium between Fe i and Fe n in r Cet. 
Yet non-LTE with pure electronic collisions cannot be excluded 
completely, because an abundance difference of A(Fe n - Fe i, 5k 
= 0) = -0.08 is only by a factor of 1.5 larger than its error bars. 

We come to the conclusion that non-LTE is a must when 
analyzing very metal-poor stellar spectra, whether hot or cool. 
LTE leads to a discrepancy in abundances between Fe n and Fe i 
of 0.09 dex for HD 84937, which is still at the level of la, and 
0.21 dex for HD 122563. As can be seen from Table[3j A(Fen - 
Fe i) should approach in non-LTE with a Sh value of between 
0.1 and 1 for HD 84937 and between and 0.1 for HD 122563. 
With S H = 0.1, A(Fen - Fei) = -0.08 dex for the first of them 



and +0.05 dex for the second one. Most F5 - K0 type stars in 
the super-solar metallicity down to [Fe/H] = -2.5 domain have 
stellar parameters T e g and log g in between those for HD 84937 
(6350/4.09) and HD 122563 (4600/1.60). Therefore, the uncer- 
tainty in the estimated Sh value is expected to result in an abun- 
dance error of no more than 0.08 dex in the non-LTE calculations 
for Fei in such a type of stars. Note that this abundance error 
decreases towards higher metallicity. For comparison, an abun- 
dance error due to using the LTE assumption can be as large as 
0.2 dex. 



6. Uncertainties in the iron non-LTE abundances 

To assess the influence of crucial atomic data and stellar param- 
eters on the accuracy of iron non-LTE abundances, test calcula- 
tions were performed for Procyon, HD 84937, and HD 122563. 
For each parameter or atomic model that we varied, we de- 
rived stellar Fe abundances from the Fe i and Fe n lines and then 
the average values. The results of the tests are summarized in 
TableU 

6. 1 . Hydrogenic collisions 

The effect of including hydrogenic collisions in our SE calcu- 
lations is shown in the first four lines of Table[4] We discussed 
above how the Fe i and Fe n abundances vary depending on the 
5h value. Here, we concentrate on the effect of using the atomic 
model, where H i collisions were taken into account for the al- 
lowed transitions, but were neglected for the forbidden transi- 
tions. For a differential analysis of stellar spectra, the calcula- 
tions were performed with Sh = 0, 0.1, and 1 not only for the 
three tested stars, but also for the Sun. As expected, neglect- 
ing hydrogenic collisions for the forbidden transitions leads to 
substantially stronger departures from LTE for Fe i in the VMP 
stars. With the modified model atom, no single Sh value lead- 
ing to consistent Fe i and Fe n abundances in both HD 84937 and 
HD 122563 was found. 

6.2. Completeness of the model atom 

For the Sun and the reference stars, calculations were also 
performed with the reduced atomic model, where all the pre- 
dicted high-excitation levels of Fe i were removed. As shown in 
Sect. 12. 31 the use of the reduced model atom leads to a stronger 
overionization of Fe i compared to that calculated with our fi- 
nal model atom. As a result, the Fe i based non-LTE abundances 
obtained with the reduced model atom are higher than the cor- 
responding values in Table [3] For the Sun, the abundance dif- 
ference amounts to 0.08 dex with Sh = 0.1. The changes in 
stellar differential abundances are smaller (TableHJi. An excep- 
tion is HD 122563, where the Fei mean abundance changes by 
+0.07 dex. We caution against applying an incomplete model 
atom of neutral iron or any other minority species to a non-LTE 
analysis of stellar spectra. 

6.3. Stellar parameters 

We investigated how the uncertainties in stellar parameters influ- 
ence our choice of the Sh scaling factor. For both VMP stars of 
our sample, the test calculations were performed with a 0. 1 dex 
(more than lcr) lower gravity and an 80 K (approximately 2<x) 
revised effective temperature, upwards for HD 84937 and down- 
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Table 4. Uncertainties in the non-LTE analysis of Fe i and Fe n in selected stars. 



Changes in non-LTE abundances, [Fe/H], relative 
to those in Col. 10 (5 H = 0.1) of Tabled] 

Input parameter HD 61421 HD 84937 HD 122563 

Fei Fen Fei Fen Fei Fen 



Hi collisions: 



Sh=0 


-0.02 


0.01 


+0.04 


+0.02 


+0.19 


+0.11 


Sh = 1 


-0.03 


0.01 


-0.13 


0.00 


-0.13 


0.00 


Sh = 2 


-0.03 


0.00 


-0.14 


0.00 


-0.15 


0.00 


Sa = 0.1/0 for 


0.00 


0.00 


+0.06 


0.00 


+0.08 


+0.03 


allowed / forbidden transitions 














Reduced model atom, 5k = 0. 1 


+0.01 


+0.01 


+0.02 


+0.01 


+0.07 


0.00 


Stellar parameters: 














Ar eff = -80 K 










-0.14 


+0.01 


Ar eff = +80 K 


+0.06 


0.00 


+0.06 


+0.03 






Alogg = -0.1 






0.00 


-0.03 


+0.02 


-0.04 


AV m ic = -0.1 kms" 1 


0.02 


0.04 


0.00 


0.01 


0.02 


0.01 



wards for HD 122563. The results of these tests can be summa- 
rized as follows. 

- The variation in T e g, log g, and V m \ c within the error bars 
does not help to achieve the Fe i/Fe n ionization equilibrium 
in HD 84937 and HD 122563 under the LTE assumption. 

- Reducing the gravity by 0. 1 dex does not change the excita- 
tion gradient of the Fe i based abundances in each of these 
two stars. A relatively steep trend of [Fe/H]i versus E exc 
obtained in HD 122563 can essentially be removed with a 
downward revision of the effective temperature by 80 K. 

- The uncertainties in stellar parameters do not influence our 
conclusion that there is a need for inelastic collisions with 
hydrogen atoms (or other thermalizing collisional processes 
not involving electrons) to establish the SE of iron in the at- 
mospheres of metal-poor stars. With the downward revised 
T e ff or log g, HD 122563 favors a low efficiency of Hi col- 
lisions with 5h = and 0.1, respectively. With the higher 
r e ff or the lower gravity of HD 84937, the difference between 
[Fe/H] i and [Fe/H]n is removed with 5h - 1. 



6.4. Stellar atmosphere and line formation modelling 

Only two non-LTE studies for iron exist in the literature that 

Shchukina & Truiillo Buenol (1200 ll) 



go b eyond the ID an a lysis. 



and IShchukina et al.l d2005l) performed 1.5D non-LTE calcu- 
lations in 3D model atmospheres of the Sun and the metal- 
poor subgiant HD 140283. They found that the departures from 
LTE are stronger in 3D than in ID model and gave a 3D 
non-LTE abundance correction of A3d,nlte = log £3d,nlte ~ 
log£3DLTE = +0.07 dex for the Su n and about +0.9 dex for 
HD 140283. IShchukinaet aljj200l have demonstrated that the 
ionization equilibrium between Fe i and Fe n in HD 140283 can- 
not be achieved with T e ff = 5700 K, log g = 3.7, and [Fe/H] 
= -2.5, and they obtained surprisingly large statistical errors 
of <Ti og£ =0.19 dex and 0.48 dex for the abundances deter- 
mined from the Fei and Fen lines, respectively. They argued 
that adopting T eS 5600 K and [Fe/H] -2.0 substantially 
reduces the discrepancies between the two ionization stages. 
We note that this effective temperature is in conflict with the 
recently im proved IRFM temperature of this star, either T e g 
= 5755 K ( Gonza lez Hernandez & Bonifacio! 120 09) or T e g = 



5777 K dCasagrande et alJl2010l) . An explanation for this result 
can be the use of an incomplete model atom of Fe i and ignor- 
ing inelasti c collisions with Hi at oms. The Fei term diagram 
applied by IShchukina et al.l (120051) is, in fact, complete up to 
E exc =5.72 eV. But, at higher energies, between E exc =5.72 eV 
an d 7.0 eV, it only c ontains about 50 % of the terms identified 
by iNave et alJ dl994h . The non-LTE+3D problem for iron still 
awaits a more satisfactory solution. 



7. Conclusions 

In this study, a comprehensive model atom for neutral and 
singly-ionized iron was built up using atomic data for the energy 
levels and transition probabilities from laboratory measurements 
and theoretical predictions. With a fairly complete model atom 
for Fei, the calculated statistical equilibrium of iron changed 
substantially by achieving close collisional coupling of the Fe i 
levels near the continuum to the ground state of Fe n. There is 
no need anymore for the enforced upper level thermalization 
procedure that was applied in the previous non-LTE a nalyses 
ilGehren et al Jl2001at BfomeT^]l2003l: ICollet et al.ll2005l) . 

Non-LTE line formation for Fe i and Fe n lines was consid- 
ered in ID model atmospheres of the Sun and five reference 
stars with reliable stellar parameters, which cover a broad range 
of effective temperatures between 4600 K and 6500 K, gravi- 
ties between log g = 1.60 and 4.53, and metallicities between 
[Fe/H] = -2.5 and +0.1. We found that the departures from LTE 
are negligible for the Fe n lines over the whole stellar parame- 
ter range considered. For Fe i, the non-LTE effects on the abun- 
dances are expected to be small for stars with solar-type metal- 
licities such as the Sun, Procyon, and /3 Vir, and for mildly metal- 
deficient stars such as t Cet: a non-LTE correction is at the level 
of +0.1 dex in non-LTE with pure electronic collisions and of a 
few hundredths, when inelastic collisions with hydrogen atoms 
are taken into account in the SE calculations with 5k > 0. 1 . 

From a differential line-by-line analysis of stellar spectra 
we found that the iron ionization equilibrium is not fulfilled in 
Procyon and /3 Vir at their fundamental parameters, 6510/3.96 
and 6060/4.11, respectively, independent of either LTE or non- 
LTE. For Procyon, an upward revision of its temperature by 80 K 
removes the obtained abundance imbalance. 
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In contrast, consistent iron abundances from both ionization 
stages were obtained in r Cet at its given stellar parameters in 
both non-LTE with Sh > 0.1 and in LTE. 

Significant departures from LTE for Fe i were found in the 
two VMP stars of our sample, HD 84937 and HD 122563. Our 
results indicate the need for a thermalizing process not involving 
electrons in their atmospheres. Since there are no accurate theo- 
retical considerations of appropriate processes for iron, we sim- 
ulate an additional source of thermalization in the atmospheres 
of cool stars by parametrized H i collisions. Close inspection of 
the ionization equilibrium between Fe i and Fe n in HD 84937 
and HD 12 2563 leads us to choose a scal ing factor of 0.1 to the 
formula of lSteenbock & Holwegerl(ll984l) for calculating hydro- 
genic collisions. The uncertainty in the estimated Sh value is ex- 
pected to result in abundance error of no more than 0.08 dex in 
the non-LTE calculations for Fe i in the F5 - K0 type stars in the 
super-solar metallicity down to [Fe/H] = -2.5 domain. However, 
the situation can be significantly worse for extremely and ultra- 
metal-poor stars. Exactly for such objects, the determination of 
the surface gravity relies in most cases on the analysis of the 
Fe i/Fe n ionization equilibrium. Theoretical studies are urgently 
needed to evaluate the cross-sections of inelastic collisions of 
Fe i with H i atoms and to search for and evaluate other types of 
thermalizing processes. 

For the Sun, the use of Sh — 0.1 leads to an average Fei 
non-LTE correction of 0.03 dex and a mean Fei based abun- 
dance of 7.56±0.09. A mean solar abundance derived from Fe n 
lines varies between 7.41+0.1 1 and 7.56±0.05 depending on the 
source of gf- values. A statistical error of 0.09 - 0.1 1 dex is un- 
comfortably high for the Sun. It is, most probably, due to the 
uncertainty in ^/-values and, in part, van der Waals damping 
constants. The problem of oscillator strengths of the Fe lines and 
their influence on t he derived solar iron abundance has been de- 
bated for decades dBlackwell et all 19951: iHolweger et al.1ll995t 
iKostik et al.ll996tlCaffau et al.l2010t) . Stellar astrophysics needs 
accurate atomic data for an extended list of the iron lines, which 
could be measured in stars of very different metallicities. We 
therefore call on laboratory atomic spectroscopists for further 
efforts to improve gf— values of the Fei and Fen lines used in 
abundance analysis. 

Using our carefully calibrated model of iron, cool stars over a 
broad range of metallicities encountered in the Galaxy can now 
be analyzed in a homogeneous way to derive their iron abun- 
dance and gravity without resorting to trigonometric parallaxes. 
The dependence of departures from LTE on stellar parameters 
and an application of the non-LTE technique to the known ultra- 
metal-poor stars ([Fe/H] < -4.5) will be presented in a forth- 
coming paper. 
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Table 5. continued. 
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Transition 


Mult 


*^exc 


logg/ 


logC 6 






logs 






^mac 


Ref. 1 


(A) 






(eV) 






LTE 


S H = 


0.1 


1 


2 


0.1 




1 


2 
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4 
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7 


8 
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10 


11 


12 


13 


5855.08 


y 3 F;-e 5 H 4 


1179 


4.61 


-1.48 
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±0.10 


7.62 

±0.11 
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±0.10 


7.54 

±0.10 


7.54 

±0.10 







4491.40 
4582.83 
4508.29 
4620.52 
6369.46 
6432.68 
5284.11 
4923.93 
5018.44 
5991.38 
5264.81 
5414.07 
5197.58 
5325.55 
5425.26 
6239.95 
6247.56 
6456.38 



b 4 F 
b 4 F 
b 4 F 



b 4 F 
a 6 S 
a 6 S 
a 6 S 
a 6 S 
a 6 S 
a 4 G 
a 4 G 
a 4 G 
a 4 G 
a 4 G 
a 4 G 
b 4 D 
b 4 D 
b 4 D 



3/2 
5/2 
3/2 



7/2 
5/2 
5/2 



5/2 
5/2 ' 
5/2 
11/2 
5/2 
7/2 
5/2 
7/2 
9/2 
1/2 
5/2 
7/2. 



7 4 F° 
z r 3/2 

7 4 F° 
*• 1 7/2 

7 4 D° 

7 4 n° 
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7 6 D° 
~ ^3/2 

7 6 D° 
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- r 5/2 



37 
37 
38 
38 
40 
40 
41 
42 
42 
46 
48 
48 
49 
49 
49 
74 
74 
74 



2.84 
2.83 
2.84 
2.82 
2.89 
2.89 
2.88 
2.88 
2.88 
3.15 
3.22 
3.21 
3.22 
3.21 
3.20 
3.89 
3.89 
3.90 



-2.76 
-3.22 
-2.34 
-3.29 
-4.25 
-3.71 
-3.13 
-1.42 
-1.23 
-3.66 
-3.13 
-3.65 
-2.34 
-3.32 
-3.38 
-3.57 
-2.43 
-2.18 



Fe ii lines 

-32.02 

-32.03 

-32.00 

-32.02 

-32.06 

-32.07 

-32.04 

-32.03 

-32.04 

-32.05 

-32.01 

-32.02 

-32.02 

-32.03 

-32.04 
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-32.00 

-32.00 
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7.48 


7.47 
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7.48 


7.48 


3.5 


RU98 


7.59 


7.59 


7.59 


7.59 


7.59 


3.7 


RU98 


7.62 


7.61 


7.62 


7.62 


7.62 


3.6 


RU98 


7.47 


7.46 


7.47 


7.47 


7.47 


3.5 


M83 


7.57 


7.56 


7.57 


7.57 


7.57 


3.5 


M83 


7.56 


7.54 


7.56 


7.56 


7.56 


3.5 


M83 


7.57 


7.57 


7.57 


7.57 


7.57 


3.5 


RU98 


7.55 


7.54 


7.55 


7.55 


7.55 


3.5 


M83 


7.56 


7.56 


7.57 


7.57 


7.57 


3.8 


RU98 


7.55 


7.53 


7.55 


7.55 


7.55 


3.4 


RU98 


7.64 


7.62 


7.64 


7.64 


7.64 


3.6 


RU98 


7.62 


7.60 


7.62 


7.62 


7.62 


3.5 


RU98 


7.59 


7.59 


7.59 


7.59 


7.59 


3.7 


RU98 


7.59 


7.57 


7.59 


7.59 


7.59 


3.4 


RU98 


7.60 


7.58 


7.60 


7.60 


7.60 


3.6 


RU98 


7.56 


7.55 


7.56 


7.56 


7.56 






±0.05 


±0.05 


±0.05 


±0.05 


±0.05 







1 References to the adopted g/- values: 
BKK9 1 = Bard et al. (1991) . BK94 = Bard & 
BPS82a = Blackwell et al. (1982a) . BPS82b = 
M83 =Moitv(1983). OWL91 = O' Brian et al. 



BIP79 = Blackwell et al. (1979) . 
Blackwell et al. (1982b) . FMW88 = Fuhr et al. 
(1991). RU98 = Raassen & Uvlings (1998) 



